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Abstract 

We discuss the recently observed excess of high lepton momentum transfer, Q 2 , 
neutral current deep inelastic scattering events at HERA in the light of supersymmetry 
with broken R-parity. We find more than one possible solution. We consider the 
possibilities for testing these hypotheses at HERA, the Tevatron and at LEP. One 
lepton-number violating operator can account for both the HERA data and the four- 
jet anomaly seen by ALEPH at LEP. 
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1 Introduction 



Recently both experiments at HERA have quoted an excess at high lepton momentum 
transfer, Q 2 , in their neutral current deep inelastic scattering (NC DIS) data [[[], 0. For 
Q 2 > 15, 000 GeV 2 , HI have found 12 events where they expect « 5 from the Standard 
Model (SM) 0. For Q 2 > 20, 000 GeV 2 ZEUS observe 5 events where they expect w 2 
from the SM In addition, HI have seen a slight excess in their charged current deep 
inelastic scattering (CC DIS) [|lj. Due to the low rate, we find it is still too early to comment 
on this. The NC DIS discrepancies could be the first hint of physics beyond the SM ||. 
There have been three main suggestions predicting an excess of events at high Q 2 at HERA: 
contact interactions @, ||, ||, leptoquarks 0, and supersymmetry with broken R-parity 

If quarks (or leptons) have a substructure this could lead to an excess of events at 
high values of Q 2 at colliders. A four-fermion operator based analysis for testing effective 
theories at e + e~-colliders was first proposed in Ref . [§J . This analysis of contact interactions 
was generalised to eP deep inelastic scattering in || where the operators are of the form 
eeqq. Contact interactions are a parametrisation of unknown physics at a higher energy 
scale (A 3> y/s) via gauge-invariant but non-renormalisable operators. At HERA, they can 
lead to an increase or a depletion of high Q 2 events, depending on the sign of the interference 
with the SM-interaction, which is a free parameter. These events should not show any peak- 
like structure, since the HERA energy is well below any supposed resonances. At HERA, 
one can attempt to extract the unknown scale of the new physics. This analysis is typically 
performed one-operator at a time. But it is also possible to consider combined effects ||. 
This should be done with the present HERA data, but we do not further consider this 
possibility here. 

Both leptoquarks and supersymmetry with broken R-parity predict s-channel resonances 
in positron-proton scattering. In Grand Unified Models, the leptoquarks are typically heavy 



and well beyond the reach of colliders [13]. However, there is a set of models which predict 
low-energy leptoquarks |14| . A systematic search for such leptoquarks was first proposed in 
Ref. 0] and subsequently much effort has been dedicated at HERA to this search ||. Within 
supersymmetry, we naturally expect all new states to have masses between 100 GeV and 
ITeV in order to maintain the solution to the gauge hierarchy problem [[n|. Depending 
on the magnitude of the Yukawa coupling and the remaining supersymmetry spectrum, 
squarks can behave just like leptoquarks at HERA. But squarks can naturally have further 
decay modes and their discussion is thus more general. In the following, we focus on the 
implications of the recent HERA data for supersymmetry with broken R-parity. When 
relevant we shall emphasise below where the leptoquark description is distinct. 

In Section ||| we briefly review supersymmetry with broken R-parity and focus on the 
relevant terms for HERA. In Section [3] we review the indirect bounds on these operators. 
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In Section ^ we determine which R-parity violating operators could possibly lead to the 
observed excess at HERA while still being consistent with the existing indirect bounds. 
In Sections |5] through |7| we discuss how these operators could be tested in future runs at 
HERA, as well as with present and upcoming data at LEP and at the Tevatron. In Section 
U we present our conclusions. In the Appendix we have collected some of the relevant 
formula for completeness. 



2 Super symmetry with Broken R-Parity at HERA 

When minimally extending the particle content of the SM to incorporate supersymmetry 
one must add an extra Higgs SU(2) L doublet and then double the particle content. The 
most general interactions of these particles consistent with supersymmetry and SU(3) C x 
SU(2)l x U(1)y gauge symmetry are those of the minimal supersymmetric standard model 
(MSSM) |16|] as well as the superpotential termsHilTH] 

KjkLiLjEk + X^LiQjDk + X'-j k UiDjDk- (1) 

Here L (Q) is the lepton (quark) doublet superfield, and D, U (E) are the down-like and up- 
like quark (lepton) singlet superfields, respectively. k = 1,2,3 are generation indices. 
The last two sets of terms in (|l|) lead to rapid proton decay The solution of this 

problem in the MSSM is to impose R-parity, 

R p = (_l)3B+L+25 ; (2) 

a multiplicative discrete symmetry. Here B denotes baryon number, L lepton numbeif] 
and S the spin of a field. All SM fields have R p = +1; their supersymmetric partners 
have R p = — 1. This solution is not unique. There are many models which protect the 
proton but allow a subset of the terms in ([I]) jlT], 20, 21 1. This subset can be as small as 



two operators even for a gauge symmetry |2U|| . All these alternative solutions are denoted 
"R-parity violation" . In the following, we shall focus on the subset of the operators ([]]) 

KjkLiQjDk- (3) 

At HERA these operators can lead to resonant squark production 

e + + Uj — > d k , (4) 
e + + d k — > Uj. (5) 

This was first proposed by J. Hewett in Ref. || where she considered the direct R-parity 

violating decay of the squark to the initial state. The processes were discussed in more 

lr rhere is the further term KiLiH u which violates lepton number. If the soft-breaking terms are universal 
it can be rotated away |l7fl . 

2 This should not lead to any confusion with the lepton doublet superfield, L, of Eq.(nh. 
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0.26^ 
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0.63^ 


0.002^ 


0.004 


0.06 


0.06 


0.13 


0.087 


0.55 


0.13 


1.3 


0.003 



Table 1: Indirect bounds on first lepton generation operators LQD. The first line is the 
bound for a scalar fermion mass of 100 GeV the second line for 210 GeV. The bounds derive 
from the following physical processes: (°) neutrinoless double beta decay p6| , ^ charged current 



universality [27, 29 1, ( c ) atomic parity violation [£g, 30, |2 



mass [31], ^ forward 

backward asymmetry [27|, and D-decays [32]. The bounds from w, ( c \ ( e \ and scale 



linearly with the squark mass (Mq/100) GeV. The bound from ^ scales with the square root 



Mg/WOGeV. The bound W scales as (M 9 /100 GeV) 2 y Mg/W0 GeV. We have conservatively 
estimated the gluino mass at M g = ITeV. We have given the 1 sigma bounds. 



detail in [|K], [11]] where the squark cascade decays via a photino were included. This enables 
the distinction from a leptoquark, depending on the supersymmetric spectrum. The full 
mixing in the neutralino sector was first considered in [[22|] for R-parity violation at HERA. 
For a scalar top quark (j = 3), the decays via the neutralino are kinematically forbidden. 
This scenario was discussed in considerable detail by T. Kon et al. |12|, |23|. The squark 
decays for the full gaugino mixing including the chargino were discussed by E. Perez [E4[ 



and in a shorter version in . This enabled a full classification of the final state topologies 
at HERA. 

Making use of all this work, we wish to determine the size and generation structure j, k 
of the coupling A' ljfc of @ which is preferred by the HERA data. We would also like to 
consider the best estimate for the squark mass. However, in order to do that we must first 
consider the indirect bounds on the couplings A' ljfc . 



3 Previous Indirect Bounds 

The operators (^|) can contribute to several processes with initial and final state SM particles 
via the exchange of virtual sleptons or squarks. Since to date all such processes have been 
observed in agreement with the SM, this leads to upper bounds on all the couplings which 
we summarise in Table |l| We give the bounds for a scalar fermion mass of 100 GeV, which is 
the standard, and 210 GeV, which is what we require below. These two bounds are related 
via scaling properties which are explicitly given in the table caption. We have updated some 



previous bounds using more recent data from the PDG []29[ . The bound from atomic parity 
violation is obtained using the "weak charge", Qw, HH an d was not previously applied 



to R-parity violation P8| . All but two of the bounds are well below the electromagnetic 
coupling. 
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There is a further set of stringent bounds from the decay K — > txvv [33]. 



\' ljk < 0.012, j = 1,2, fc = l,2,3; M = 100 GeV, (6) 
A' 132 <0.19, M = 100 GeV. (7) 

These limits scale linearly with the squark mass. However, these bounds are model depen- 
dent. If the absolute mixing of the quarks (not squarks) is purely in the down-quark sector 
then these bounds apply. If it is purely in the up-quark sector then these bounds revert to 
those of Table 1. In order to keep an open mind to all possible solutions we do not further 
consider these model dependent bounds. 

If there is a charged, first or second generation doublet slepton (e^, lighter than the 
top quark there is an additional bound from top quark decays |33| 

A' 132 < 0.4, M| = 100 GeV, (8) 

which is more strict than the one given in Table 1. 



HI have performed a direct search for supersymmetry with broken R-parity pi |. The 
bounds from this search depend on the neutralino mass. For a neutralino mass M^o = 
40 GeV < Mg they obtain 

A' yi <0.2, Mg = 200 GeV, M^o = 40, (9) 
A' 132 <0.22, Mg- = 150 GeV, M^o = 80. (10) 

In particular the last bound is significantly better than the indirect bound of Table 1, but 
is also model dependent. 

4 HERA's high Q 2 excess interpreted as s-channel squark 
production 

4.1 R-parity Violating Squark Decays 

We would now like to interpret the observed excess at HERA in terms of supersymmetry 
with R-parity violation. For this, we combine the HI and ZEUS data in order to compare it 
more easily with different R-parity violating models. For Q 2 > 20, 000 GeV, HI and ZEUS 
see a total of 10 events Jl|, 0, where 4.08 ± 0.36 events are expected from SM contributions 



35|. The total integrated luminosity of the two experiments is 34.3p6 1 , which translates 



into an excess cross-section over the SM expectation of 

a ex (Q 2 > 20, 000 Gel/) = (0.17 ± 0.07 (stat)) pb. (11) 
In Figures |l] and § we show this combined excess cross section as a horizontal band. 
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The combined reconstructed values of the hypothetical squark mass M q (where Mq is 
related to Bjorken-x by M| = xs and s is the centre-of-mass energy squared) show some 
spread between the two experiments (M q = 200 GeV for HI, M- q = 220 GeV for ZEUS). 
Combining the two experiments, we obtain^ 

Mq = (210 ± 20) GeV, (12) 

as our best estimate. 

In order to determine the contribution from R-parity violation we consider one of the 
operators at a time while assuming the others are negligible. For a given non-vanishing 
operator X'^LxQjD^ the produced Uj, or dk squark can have many decay modes, depending 



on the mass spectrum of the supersymmetric particles [10, 11, 24, E5L pqj. For all SUSY 



spectra, the squark can decay via the operator itself resulting in the interactions 

e + + Uj — > dk — > e + + Uj, (13) 
e + + Uj — > d k — > z/ + c/j, (14) 
e + + — > — >■ e + + c? fc , (15) 

at HERA. The first and third have equivalent initial and final states to NC DIS and can 
thus contribute to the observed excess. Analogously, the second can contribute to CC DIS. 
Let us for now assume these are the only decay modes of the produced Uj or dk squark. We 
then compute the production cross section using |9|, [iT| and the MRSG structure functions 



37 1 . For (|13|) and (14) we include the extra contribution to the width for the other decay 



mode. The correction to the narrow width approximation by using the full resonance width 



is ~ 10%. 



In Figure [I], we plot the production cross section <r(e + + dk — > Uj — > e + + dk) for 
Q 2 > 20,000(^6^ as a function of the R-parity violating couplings A' ljX , A' lj2 , A' lj3 . For 
each coupling, the three curves are for the squark masses M^. = 200, 210, and 220 GeV 
from top to bottom, respectively. The cross section is largest for X' l j l due to the incoming 
valence quark. The other two cross sections are suppressed because of the significantly 
smaller strange and bottom sea-quark structure functions s(x), b(x). 

In Figure §, we analogously plot the cross sections <r(e + + Uj — » d k — » e + + Uj) as a 
function of the couplings X' nk and A' m for the three squark mass values Mj fc = 200, 210, 
and 220 GeV. Here, we have included the full R-parity violating width from the two dj 
decay modes. These cross sections are both suppressed because of the small up and charm 
sea-structure functions u(x) and c(x), respectively, as well as the increased width. There 

3 The errors quoted on the mass of the excess events (M e ) in Table 6 Q is not the actual error on the 
reconstructed mass of any supposed resonance. This error is most likely larger as can be seen from the 
broad width of the distributions in Fig. 2b of Ref. Analogous, remarks apply to the ZEUS analysis 0. 
We therefore feel justified in combining the data in this way. 
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is no contribution for A' 13fe because of the unknown and suppressed top-quark structure 
function, t(x). 

We now compare the plotted cross sections with the hatched band of Eq. (p]) , and the 
bounds in the second row of Table |l[ We obtain the following set of solutions where R-parity 
violation can explain the excess in the HERA data 



(16) 



In each row we first present the approximate value of the required Yukawa coupling, then 
denote the produced squark and the mass range of the squark which is viable. Thus we 
obtain many solutions; however, the preferred solutions which are well within the constraints 
are ([16|.l), (|16|.4) and ([16|.5) with a produced scalar charm quark and a scalar top quark 
("stop"), respectively. Before discussing the other squark decay modes, we note that within 
supersymmetric unification scenarios there is a possibility for a very light stop, possibly even 
the LSP fl38| . For a stop-LSP, the process (|1^) would be the only decay. 
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There is one possible further solution, which is intriguing. Combining the solutions 
flTB|.2) and (|TB|.3), it is possible that HERA has produced two different squarks: a scalar 
charm and a scalar bottom. The required coupling is then reduced by a factor \/2 to 
A' 123 ~ 0.28, which is well away from the indirect bound. 

(7) A' 123 = 0.28, c & b, Mg « 210 GeV (17) 



We make one further point []39 |. The leading order QCD correction to the resonant 



squark production is presently not known. On-shell squark production is effectively a 
2 — > 1 process for which the QCD corrections are typically positive and large ||y]. If this is 
indeed the case for resonant squark production at HERA this would give additional lee- way 
in the Yukawa-coupling for the more marginal solutions. The required coupling would be 
reduced to A'/ \f~K, where K is the QCD if-factor. 



4.2 Gaugino Decays of the Squark 

In this Subsection we wish to discuss the case where the produced squark is not the LSP 
and can thus have further supersymmetric decay modes. (This section has no analogy for 
leptoquarks.) As noted above, supersymmetry dramatically extends the spectrum of the 
SM. It is way beyond the scope of this paper to perform a systematic study of all the 
possible decay chains. Instead, we focus on two specific cases which are well motivated by 
the renormalisation group studies of the MSSM |4l| . (a) The lightest neutralino, \ii is 
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Figure 1: Excess cross-section a(e + dk — ► e + dk) for a s-channel Uj-squark resonance for Q 2 > 
20, 000 GeF as a function of the coupling X' ljk for the three values of M fij = 200, 210, 220 GeF 
(top to bottom). The solid lines show the cross-sections for a non-zero coupling A^ (i.e. the 
d valence quark contribution), while the dashed-dotted and dotted lines show the cross-sections 
for the non-zero couplings A^, and A'^ (i.e. the s, and b sea quark contributions), respectively. 
Here, we assume that the only allowed squark decay mode is Uj — ► e + + dk- The hatched region 
shows the high Q 2 excess cross-section a ex = (0.17 ± 0.07) pb. 
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Figure 2: Excess cross-section a(e + Uj — > e + Uj) = a(e + Uj — > vdj) for a s-channel d^-squark 
resonance for Q 2 > 20, 000 GeV as a function of the coupling A' ljfc for the three values of = 
200, 210, 220 GeV (top to bottom). The solid lines show the cross-sections for a non-zero coupling 
A' llfc (i.e. the u-bar sea quark contribution), while the dotted lines show the cross-sections for the 
non-zero coupling A' m (i.e. the c-bar sea quark contributions), respectively. Here, we assume that 
the only allowed squark decay modes are — > e + + Uj and dj~ — > v + dj. The hatched region 
shows the high Q 2 excess cross-section a ex = (0.17 ± 0.07) pb. 
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Figure 3: Cross-section a{e + d — > x°c) f° r single neutralino production at HERA for the coupling 
X'i2i =0.1 and Mg = 210 GeV as a function of the neutralino mass. The solid, dashed and dotted 
lines show the cross-sections for a neutralino which couples purely Wino-, Bino-, or Higgsino-like, 
respectively. 
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Figure 4: Effect of the additional decay mode c — ► x°c on the excess cross-section o ex (e + d — » 
e + d) for a s-channel c-squark resonance for Q 2 > 20,000Gey, A' 12 i = 0.1 and Mg = 2WGeV, 
as a function of the neutralino mass. The hatched region shows the high Q 2 excess cross-section 
a ex = (0.17 ± 0.07) pb. 
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Figure 5: Cross-section a(e + d — ► for single chargino production at HERA for the coupling 
A' 12 i = 0.1 and = 210 GeV as a function of the chargino mass. The solid and dashed lines 
show the cross-sections for a chargino which couples purely Wino-, or Higgsino-like, respectively. 
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lighter than the produced squark, M^o < M^, M^ k and (b) The lightest chargino, xt ■> is 
lighter than the produced squark, M~+ < M^, M^. 
(a) If 

M^< M„ r M, lr (18) 
then we have the following additional interactions at HERA 

e + + % ->J fc + (19) 

e + + 4 + (20) 

In the second process, the case j = 3 is suppressed by the large top-quark mass and we do 
not further consider it. 

The production cross section depends on the admixture of the neutralino mass eigen- 
state. In order to determine the cross section we focus on three special limiting cases! a 
purely Wino-, a purely Bino- and a purely Higgsino-x5- I n Figure^ we plot the production 
cross section a(e + + dk — > Xi + u j) as a function of the neutralino mass for these 3 special 
casesfj We have fixed the squark mass at 210 GeV and the Yukawa coupling X' 1 j 1 = 0.1. 
For a purely Higgsino-neutralino, the production cross section is very small (< 10~ 6 pb) and 
is not observable at HERA. For a gaugino-neutralino, the cross section is ~ lpb for neu- 
tralino masses as large as 150 GeV or 180 GeV. With the present luminosity, HERA could 
possibly have already produced several tens of neutralinos, upto about 30 per experiment. 
The neutralino decay modes are |p22 



Xi — > {e Ujd k , e + ujd k , v e djd kl v e djd k } ■ (21) 

The branching fractions are not simply 0.25 each for massless final states but sensitively 
depend on the admixture of the neutralino [22 1. The most visible decay for a positron beam 



is x\ — > e~Ujdk, i.e. to the final state electron [II|. This requires charge identification of 
high pt electrons and can be searched for in the present and upcoming data. 

By including the neutralino decay of the squark we have increased the total decay 
width of the squark. This in turn reduces the resonant production cross sections plotted 
in Figure [l|. In Figure |], we show this effect for the three neutralino admixtures and for 
= 0.1, M(v,j) = 210 GeV. The change in the leptoquark-like cross section (squark 
LSP) is negligible for the Higgsino-neutralino. For a Bino-neutralino, the cross section can 
drop by about a factor two. For a Wino-neutralino the production cross section can drop 
by more than an order of magnitude. 

If we include a Bino- or Wino-x? which is lighter than the produced squark, we must 
reconsider the solutions presented in Section 4.1. For the same coupling A' ljfc , the rate for 
the reactions (|T^,|I^) will go down when including the neutralino. This can be compensated 
by increasing the Yukawa coupling, if allowed. It is only for the solutions (|16|.lJT^.4,p^.5) 

4 The single neutralino/chargino production cross-sections are quoted in the Appendix for completeness. 
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and (17) that this is possible. For example, for solution (|l6|.l), producing a scalar-charm 
quark, we can raise A' 121 = 0.04 by a factor 2.5 to A' 121 = 0.1. We then obtain a new solution 
for a purely Wino-x? with mass M^o = 125 — 165 GeV. This can be seen in Figure |j. The 
branching ratio for the direct decay c — > e + d is reduced to about 0.1. For a purely Bino- 
x\ only a modest increase in coupling can be allowed, otherwise the suppression is not 
sufficient. 

We can similarly consider the effects of the decay d k — > x\ + dk on the production cross 
section a(e + Uj — > e + Uj). The SU{2) singelt (i-squark only couples significantly to a purely 
Bino-Xi- This only leads to marginal affects, as seen above. Except for the two squark 
solution (0), these solutions are already marginal and we do not further consider them. 

(b) We can repeat the analogous analysis for a chargino with 

M-+<M Sj ,M 4 , (22) 
in which case we obtain the additional interactions 



e 



+ 



+ Uj -> d k -> Xi + u k, (23) 
e + + d k ^uj ^xt + dj. (24) 

Here we also allow for the case j — 3. We consider the two limiting cases, where the chargino 
is pure Wino-j^ or pure Bino-x^- In Figure [5] we plot the single chargino production cross 
section a(e + d k — > xfdj) for these two cases for X[j 1 = 0.1, and M{iij) = 210 GeV. For 
the pure Higgsino case, the cross section is again highly suppressed. These decays can be 
neglected. For the purely Wino case, a(e + dk — > xfdj) ~ Ipb, for M(xt) < 180 GeV. 
This suppresses the NC-like production as in the neutralino case and we obtain additional 
solutions. This is particularly relevant for t production. 
The chargino decays to 

xt -> {is e Ujd k , u e Ujd k , e~djd k , e + djd k } ■ (25) 

Again, the best search mode for a positron beam is the final state electron which can be 
detected at HERA. 



Summarising, the solutions (|TJ.l,|Tg.4,|T6[5) as well as the special solution (\L7\) can allow 
for substantial decays to gauginos. The other solutions can only allow for Higgsino-like 
gauginos to be lighter since the decay rates are negligible. Note that the new processes via 
gaugino decays to positrons and neutrinos would be reconstructed at random x and Q 2 for 
both NC and CC. 

5 Tests at HERA 

There are several tests of the R-parity violating hypothesis which can be performed at 
HERA. 
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1. If HERA is switched back to e p collisions the processes flUf) -(0) change to their 
charge conjugate 



e~ + Uj — > d k — > e _ + itj, (26) 
e~ + Mj — > dk^Ve + dj, (27) 
e~ + rf fe — ► Uj — > e~ + d k . (28) 



For = 1 or j — 1 the event rates would increase or decrease, respectively by |37 
u(x,Q 2 )/u(x,Q 2 ) ~ 100, 

for x = 0.45, Q 2 = 4- 10 4 GeV 2 . (29) 

d(x,Q 2 )/d(x,Q 2 ) « 0.03, 

Using the e~P data from earlier runs at HERA, the first case would exclude any 
solution with an incoming w-quark. But no such solution was found due to stringent 
low-energy bounds. For the solutions (|TB[1,|T^[4) this change should lead to a non- 
observation. For the other solutions with higher generation incoming quarks, there is 
no difference between the q(x, Q 2 ) distribution and the q(x, Q 2 ) distribution and we 
expect no effect. 

2. In the operator LiQjD k , L\ refers to the first generation lepton SU(2) doublet su- 
perfield. Therefore, the positron in (|13D-(|15|) is right-handed, and the electron in 
(|2T)|)-(|2"8"|) is left-handed. Thus, when the electron/positron polarisers are installed the 
event rate should double or vanish depending on which polarisation is chosen for the 
lepton beam 0. This effect is independent of the quark generation indices. 

3. Within supersymmetry with broken R-parity, the nature of the lightest supersymmet- 
ric particle (LSP) is unknown. The decay spectrum of the produced squarks depends 
on the supersymmetry spectrum as a whole. In the previous Section, we considered 
the case where a neutralino or chargino is lighter than the squark. This leads to 
an additional electron (positron) signal for a positron (electron) beam, which can be 
searched for. 

The most important conclusion from this is that HERA itself can determine the nature 
of the observed effect. The second conclusion we draw is that it is essential for HERA to 
first continue running in the present mode and not switch to an electron beam or run at 
lower energy. This way both experiments can establish whether there is a genuine effect or 
whether it is merely a statistical fluctuation. If HERA did switch to an electron beam and 
the effect is due to the production of up-like squarks from down quarks the excess would 
be decreased and require longer running time to establish its nature. We could then end 
up in one year with two separate non-significant excesses and being no-more the wiser. 
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6 Signals at LEP 2 



In this Section we discuss three tests of the R-parity violation hypothesis at LEP 2. (1) Pair 
production of selectrons. The pair production of squarks at LEP 2 with rriq = 200 GeV is 
kinematically prohibited. Selectrons on the other hand, could be kinematically accessible 
at LEP without having been seen at HERA. In particular, they can decay through the same 
R-parity violating operator LiQjD^. With the present data samples of L ~ 25 pb~ x per 
experiment at y/s = 130 — 172 GeV, selectron masses rug < 70 GeV are accessible ||42|| . (2) 
Gauginos, if light, can be pair produced and subsequently can decay via LiQjDk to visible 
leptonic final states. (3) Virtual t-channel exchange of squarks would give a contribution 



to the SM process e + e — > qq [43, 44|. We investigate these three effects in the following. 



6.1 Selectron pair production and the ALEPH 4-jet Anomaly 

If the selectrons are lighter than the gauginos, the dominant decay channel for the Yukawa 
coupling size of interest is: 

ei,R -> Ujd k . (30) 

Selectron pair production would then give rise to 4-jet signals at LEP 2 through the process 
&l&r {ujdk){ujdk)- This scenario has been investigated in [[RJ, where it was found that 
selectron masses of = 58 GeV, m e - B = 48 GeV could explain the anomalous invariant 



mass peak of 4-jet final states observed by ALEPH ||46 |, as well as the apparent difference in 
mass AM ~ 10 GeV. It must be pointed out that the ALEPH anomaly is not seen by the 
other three LEP experiments, and its origin is at present not understood. In this Subsection, 
we will assume that the ALEPH effect is due to new physics. It is then intriguing to ask 
whether the HERA excess when interpreted as s-channel squark production is compatible 
with the above interpretation of the ALEPH 4-jet signal. In order to explain the ALEPH 
data there are several requirements on A' ljfc |45| . 



(a) The solution to the ALEPH 4-jet anomaly requires the associated production of an 
SU(2) singlet selectron (£r) and an SU(2) doublet selectron (e L ) with different mass. The 
en can not decay directly via the operator LiQjDk, only via the mixing with the e^. If this 
mixing is small then for small A' the cr will be long-lived. This in turn can be observed in 
the ALEPH detector but wasn't. Thus the coupling must be large: A' ljfe > 0.01 |45| . 



(b) The ALEPH 4-jet data excludes final states which contain b-quarks, and favours 
u,d,s quark final states; charm quarks are disfavoured, but not excluded. Including the 
tight constraints on A' m one arrives at three solutions: A' 112 , A' m A' 122 (c./. Table 0). The 
solution with A^ 2 ~ 

0.01 is favoured by the ALEPH data g|]. 

As discussed in Section 4, the HERA high Q 2 excess, too, can be explained by one of 
the three couplings (c.f. Table |3|). 

In order to fit the ALEPH data comfortably, the coupling A' 121 should be as large as 
possible. This evades the effects of lifetime (which are now amplified by the presence of 
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Coupling 


Process which can explain 
the ALEPH 4-jet anomaly 


A' 112 ~ 0.01 
A' 121 > 0.01 
A' 122 > 0.01 


— > uuss 
eL&R ^> ccdd 

&L&R CCSS 



Table 2: Allowed couplings which can explain the ALEPH 4-jets. 



Coupling 


Process explaining the 
HERA high Q 2 anomaly 


Additional signals 
at HERA 


A' m > 0.04 


e + d — > c — > e + <i 


e + <i — > c — > ( or 



Table 3: Allowed coupling which can explain the HERA high Q 2 anomaly. 



charm in the 4-jet final states). But for a large coupling (say A' 121 = 0.1), the HERA 
excess cross-section of 0.17 pb only fits the predicted cross-section if the squark coupling 
to the neutralino is large. The additional decay mode c — > x°c must reduce the predicted 
cross-section a(e + u — > c — > e + w) from 0.93 pfe to 0.17p6 (Figure f§). A neutralino with 
Mco = 125 — 165 GeV which couples with a dominant bino component indeed fits the data. 
The additional process e + u — > c — >• x°c has a cross-section of 1.4 p6 (Figure |^), which would 
predict around 50 events between the two experiment HI and ZEUS in the 1995/1996 
data; a signal which should be observable. That cross-section decreases, as the value of A' 121 
is decreased, since the neutralino mass has to be increased correspondingly to fit the dataf]. 

We conclude that the ALEPH anomaly and the HERA excess can in principle be simul- 
taneously explained by the coupling A' 121 > 0.04, although this scenario is less favoured by 
the ALEPH data. 

6.2 Other Signals at LEP 2 

Let us now discuss the constraints placed on the supersymmetric spectrum by the solutions 
of Eq. (|1~6"D . We find four scenarios which determine the gaugino spectrum: 

(a) Solution (|i"6|.l) can accommodate gauginos which are substantially lighter than Mg 
and couple electroweakly to the squarks (e.g. non-Higgsino like). 

(b) Solution (0.4) requires M^o > M~ t - M t , and - .\/, - M,,. 

5 It was pointed out in [ [si] that cancellations in the gaugino coupling of the neutralino to the left-handed 
charm squark can occur in large regions of Supersymmetric parameter space. In this case decays of the 
squark to the chargino - which do not have the same cancellations in coupling - would play the same role 
as the decays to the neutralino discussed above. 
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(c) All other cases require either large gaugino masses (~ 200 GeV), which is of no interest 
to LEP 2, or 

(d) they can accommodate light gauginos which a small coupling (e.g. Higgsino-like) to 
the squarks. 

In summary, we conclude that cases (a), (b) and (d) are interesting for LEP 2 since charginos 
and neutralinos could be discovered up to the kinematic limit of yfsjl. LEP 2 itself however 
has no way of testing the R-parity hypothesis of the high Q 2 HERA excess, since the non- 
observation of R-parity violating SUSY at LEP 2 cannot rule out that the HERA effect is 
indeed a sign of a s-channel squark resonance. 



6.3 Virtual Effects at LEP 2 from t-channel Squark Exchange 

We quote the cross-section for this effect, which has also been studied in the literature 



|3L Hi | , in the Appendix. The magnitude of the effect is small, since the main contribution 
to the cross-section is proportional to A' 2 . Figure [| shows the contribution of t-channel 
c and s exchange for Mq = 200 GeV at y/s = 190 GeV including Initial State Radiation 
corrections. Note, that the effect due to a c exchange gives a positive contribution to 
(j(e + e~ — > dd), while the s exchange gives a negative contribution to a(e + e~ — > uu). 

The overall effect on the total cross-section is shown in Table f|. We have also included 



a column with cross-sections including an anti-ISR cut Js'/s > 0.9, where s is the center 
of mass energy of the incoming lepton beams, and s' is the center of mass energy of the 
outgoing qq pair. The cut enhances the contribution of the t-channel squark exchange to the 
total qq cross-section. Nevertheless, the effect is still less than 3% for a coupling A' = 0.2. 
LEP 2 can therefore only probe large Yukawa couplings. For an integrated luminosity of 
400p6 _1 (expected to be delivered to the four experiments for the 1997 run) LEP 2 would 
see a three sigma effect on the total qq cross-section for couplings larger or equal than 
A' ~ 0.2(0.4) for s(c) squark exchange respectively. 

Very recently, the OPAL collaboration have performed the measurement |47| . For squark 
masses of 200 GeV, they are sensitive to couplings A' ~ 0.4, which is already on the verge 
of excluding the large coupling solutions in (|i~6|). 



7 Signals at the Tevatron 

There are several potential tests of the above 'supersymmetry with broken R-parity' hypoth- 
esis (produce one squark via operator LxQjD}.) at the Tevatron: (a) t-channel squark ex- 
change interfering with Drell-Yan production, (b) single squark production, and (c) squark 
pair production. 
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V irtual effect on a(eV^qq) at LEP2(Vs= 1 90GeV ) 




J I L 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos('tf) of quark q 



Figure 6: Differential cross-sections da/dcosO for M 5 ,M S = 200 GeV and A' = 0.2 respectively. 
We have included Initial State Radiative corrections, and used the anti-ISR cut y/s'/s > 0.9 (see 
text). 



19 



Anti riit 


(J I C O * CfctX 1 


(J l c c r 1 


H TT/^r i 't" nn fnfa 1 pthcc contin n 




SM nnlv 

kjiv± win y 


SlA/T Flint; f-PiiP'riTipl pvpn pi ti g"p 
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u ye e — * (/(/; 




18 62Db 


18 73Db 


4-0 12% 


*/s'/s > 9 


3 31nb 


3.42pb 


+0 53% 




<r(e + e~ — > uu) 


<r(e + e~ — > uu) 






SM only 


SM plus t-channel 5 exchange 






18.40pb 


17.77pb 


-0.68% 


/s'/s > 0.9 


5.30pb 


4.77pb 


-2.58% 



Table 4: Cross-section values (including Initial State Radiation corrections) for the SM process 
and the SM plus t-channel squark exchange for Mq = 200 GeV and A = 0.2 at = 190 GeV. 



(a) The operators LiQjD can also contribute to Drell-Yan production via the t-channel 
exchange of a squark. The dominant effect is the interference with the SM. This has been 
studied in pEBfl . However, the above study did not consider the for us relevant range at the 



Tevatron, basically because it was not feasible with the projected luminosity. In light of 
the planned Tevatron upgrade, we have repeated the analysis for our solutions. However, 
even for the large couplings we do not find a measurable effect (< 5%) in this case, 
(b) The scalar quarks can be singly produced via the parton-level processes 

g + uj -> e + + 4, (31) 
g + d k — > e~ + Uj, (32) 

as well as the complex conjugate production mechanisms. Here g denotes a parton gluon. 
This is completely equivalent to single leptoquark production at the Tevatron which has 
been analysed in f4~8fl . The production cross section is proportional to \'^ k . For \\- k = 0.1 
and Mlq = 200 GeV, we have a = 0.02 pb' 1 ji8] for an incoming first generation quark. 



Thus for our preferred solutions with A' 121 , A' 131 = 0.04 (0.1) we do not expect any observable 
signal. For our more marginal solutions we have incoming sea-quarks which significantly 
suppresses the cross section. Thus even for A' = 0.3, 0.4, we do not expect an observable 
cross section. 

(c) Squarks can be pair produced at the Tevatron via0 

g + g -> ? + (33) 

q + q — > q + q. (34) 



6 The second process includes t-channel electron exchange via the R-parity violating operator. For the 
couplings we consider this contribution is negligible. 
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If the squark only decays via the dominant first lepton generation R-parity violating opera- 
tor LiQjD k then this is equivalent to first generation leptoquark production and decay. DO 
have recently updated their analysis [|9| and obtained preliminary lower bounds on such a 
leptoquarkQ 

M LQ > 175 GeV, for BR{(f> LQ e±g) = 1, 
DO: M LQ > 147 GeV, for BR((j) LQ -+ e±g) = 0.5, (35) 
M LQ > 71 Ge^, for 5^(0^ -> e±g) = 0.0. 

The first bound applies for a two-body leptonic branching ratio of 1, i.e. the only decay 
mode is 0lq — > e 1 * 1 + q. This can immediately be reinterpreted as a lower bound for up- 
like squarks, provided they can not decay via gauginos or other supersymmetric particles. 
The second, weaker bound applies to the case where there is a second decay without a 
charged lepton and which has equal decay width to the charged lepton decay. This can be 
reinterpreted as a lower bound for d-like squarks, since as shown in Eqs.(|l^,0) they have 
two equal decay modes (provided the final state quark masses can be neglected). 

Before applying these bounds we note that they only employ the tree-level cross-section 
for leptoquark pair-production. In our case, for squark pair production, the next-to-leading 



order calculation has been done. The production cross section is increased by « 2 ||50|| . 
If we include this in the DO-plot the more stringent bound increases by about 20 GeV to 
195 GeV. This should be further strengthened by combining the DO and CDF leptoquark 
bounds, and can most likely significantly cut into our required mass range (|I2"D . With the 
Tevatron upgrade the entire mass range can be covered. 

For the down-like squarks, which have a decay branching ratio to charged leptons < 0.5 
these DO bounds are not relevant and our previous analysis of Section 4.1 applies. For the 



up-like squarks with only R-parity violating decays, solutions ([16| . 1 JTq. 2 JTq. 4,[L6|. 5) and (|Ti 
these bounds apply directly and are relevant. At present they do not exclude any solution, 
but after combining the numbers of DO and CDF the low-mass range of our solutions could 
be excluded. In this case one must include the gaugino decays of the squarks. As discussed 
in Section 4.2, the branching ratio BR(c — > e + + d) can be reduced to 0.1 for A' 121 = 0.1 and 
a neutralino mass ~ 165 GeV. In this case there are no relevant bounds from the Tevatron. 



8 Conclusions 

We have discussed the recently observed excess at HERA in the high Q 2 data in the light of 
supersymmetry with broken R-parity. We have two solutions which are in good agreement 
with all the data: 

• ^i2i > O-O^ and the production of a scalar charm quark, c. If the c is the LSP this 
solution can most likely be tested at the Tevatron by combining all the existing data, 
7 These numbers are from March, 1997. 
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otherwise with the upgrade. If it is not the LSP, then the c can decay to a neutralino 
or a chargino. The resulting spectrum can not be tested at the Tevatron with the 
leptoquark pair production search. The neutralino production however should be 
testable at HERA with a cross section of 0(lpb). The signal is a charged lepton of 
opposite sign to the beam lepton. 

• A' 131 > 0.04 and the production of a scalar top quark, stop, t. The possible solutions 
are analogous to the scalar charm with the exception that the decay to the neutralino 
is kinematically prohibited in most cases. In the renormalisation group studies in the 
MSSM a chargino LSP is disfavoured. Again the stop-LSP solution can be tested at 
the Tevatron and for a light chargino at HERA via chargino decay to an electron or 
positron. 

We also have a set of further solutions which involve larger values of the coupling. If 
some of the model dependent bounds we have discussed indeed hold, these solutions are 
excluded. For the model independent indirect lower bounds these solutions are still viable. 

• ^123 = O-^ an d t ne production of scalar charm quarks, c. This solution prefers the 
lower mass range of the resonance solution to the HERA data (< 210 GeV). If the 
c is the LSP, then the dilepton search at the Tevatron can possibly test this model 
with the present data, if both experiments are combined. It can most certainly test 
it with the upgrade. 

• ^'123 = O-^ an d the production of scalar bottom quarks, b. This solution also prefers 
the lower mass range (< 210 GeV). It can not be tested at the Tevatron since the 
charged leptonic branching ratio is less than 1/2. However this solution can be tested 
by the t-channel process at LEP 2. 

• A' 132 = 0.3 producing a t squark. This solution allows for the full mass range. As 
for the c, this model can be tested at the Tevatron if the squark is the LSP. If it is 
not the LSP there is room in the coupling to allow for stop decays via the chargino 
(neutralino decays are kinematically forbidden or suppressed by the top quark). 

We have one further solution: A' 123 = 0.28 and the production of two kinds of squarks, 
a c and a b. The tests are analogous to those discussed above for the separate solutions 
(A' 123 , c) and (A' 123 , b) . 



It is amusing to note that in Ref.pO| models with an anomaly-free family dependent 
U(1)r i?-gauge symmetry were constructed which predicted the dominant R-parity violating 
operators L 1 Q 1 Z) 2 , and L 1 Q 2 Z) 1 . 

We eagerly await further data to confirm or reject this hypothesis. It might then also 
be possible to make a statement about CC DIS and contributions from supersymmetry 
decays. 
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9 Appendix 

9.1 Single Neutralino Production Cross-section at HERA 



The cross-sections for the process e P — ► }^ + X in the approximation where the neutralino 



is a pure photino have been previously calculated in |TT[ . In this Section, we generalise the 
results to the case of a neutralino (which is an admixture of the photino, wino and higgsino 
weak eigenstates). The differential cross-section for single neutralino production via the 
R-parity violating coupling A' ljfe at HERA may then be written as 

d (x,Q 2 )\M(e + d k ^ x \)\ 2 ^ 



dxdQ 2 16irx 2 s 2 

q u (x,Q 2 )\M(e + u J -, X °d k )\ 2 ) (36) 

where qd(x,Q 2 ) and q u (x,Q 2 ) give the probability of finding a o^-quark and t^-quark re- 
spectively inside the proton; s, x, Q 2 are the center of mass energy, the Bjorken scaling 
variable and the momentum transfer squared. The Matrix elements can be obtained from 



22 1 by crossing and are given by (neglecting initial state masses): 



\M{e+d k ^ X \)\ 2 = ^{ 



-—————[(a( Uj f + b{u) 2 ){s - m 2 x o - m 2 uj ) 

'"'uj) ~r L uj UL uj 



+4a(uj)b(u)m U jm x o] 
(u - mi, J 

v dk' 

s ~ mlj 

[(a - mlj) 2 + r^-m^][t - m\ 



■[2a(uj)b(e)m U jm x os 
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+b(e)b(u){s 2 + t 2 -u 2 + {m 2 x o + m 2 uj )(u - s - t))\ 
(t 2 + u 2 -s 2 + (m 2 x0 + m 2 uj )(s -t-u)+ 2m^mJ )] 

2 

g — ffi't . _ 

[2a(uj)b(d)m U jm x os 



[(* - m lj) 2 + r lj m lj] [u - m 2 & ] 
+b(u)b(d)(s 2 + u 2 -t 2 + (m 2 x0 + m 2 UJ )(t - s - «))]} (37) 

where s, t, u are the Mandelstamm variables defined as s = (p e +Pdk) 2 = %s, t = (p e —Px ) 2 = 
—Q 2 , u = m 2 + m 2 dk — s — t; g = sin e g ; m x o, m U j are the masses of the final state particles, 
and me,maj,m^ k are the masses of the exchanged scalar SUSY particles; Taj is the total 
width of the scalar uj. 

|Me + %-X°4)| 2 = ^{ 



(u - mlj) 

-4a(d fc )6(J)m dfe m x o] 
m 2 , — t 



;b(e) 2 (m 2 -t) 



it - ml) 2 

— 7T-. 7T-b(e)b(u)\t 2 - s 2 + u 2 + (m 2 + m 2 dk )(s - t - u) + 2m 

(t-m|)(w — to|) x 

— dfc [2a(d k )b(u)m dk m x vs 

-b(u)b(d)(u 2 + s 2 - t 2 + (m 2 + m^)(t - s - it))] 

[( s - m | fc )2 + r| fc m| fc ][*-m|] 

-6(e)6(rf)(s 2 + t 2 - u 2 + (m 2 = m 2 dk )(u - t - s))]} 



where s,t,u are now defined as s = (p e + p U j) = xs, t = (p e — p x o) = —Q , u 
m 2 o + mlj — s — 
the x, Q 2 range 



m 2 o + m 2 • — s — t. The total cross-section can be obtained from Eq. fl3~6|) by integrating over 



Xmin = (m 2 +m 2 J/s 

•Emax 1 

Qmin = S - m 2 

Qmax = * (39) 

and m/ s is the mass of the final state quark, i.e. m U j, rridk for Eqns. (|37|) ,(|38|) respectively. 
The coupling constants a,b are given in Table A.l of reference EI. The cross-section for 
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the process e P — > x° + X can be obtained from the above result by charge conjugation: 

q u (x,Q 2 )\M(e- Uj ;^X°4)| 2 ) (40) 
where \M(e~d k — > x°%)| 2 = |.M(e + <4 — > x ^)! 2 and |.M(e~'Uj — > X°c4)| 2 = |A^(e + % — > 

x°4)| 2 - 

9.2 Single Chargino Production Cross-section at HERA 

The differential cross-section for single chargino production via the R-parity violating cou- 
pling X[ j k at HERA may be written as 

?„(x,Q 2 )|JW(e+u J ^X + «*)| 2 ) (41) 
where the Matrix elements can be obtained from by crossing and are given by 

\M{e+d k ^ X +d 3 )\ 2 = ^{ 

^;{ll + ll){ml + -t) 



\ b ln uj) i 1 uj' n uj 



* - ;Ke{ lL 8* L [s 2 + t 2 -u 2 



[(a - mfy) 2 + r^m^.] [t - m?] 
+(m| + m 2 + )(w - s - t)] + 2 lL 8* R m dj M x+ s}} (42) 

where s, t, u are defined by s = (p e +Pdk) 2 = xs, t = (p e —p x +) 2 = ~ Q 2 ■> u = rn^+m^ — s—t. 
And 



g 2 \' 2 - 1 

\M{e + u 3 -> x + w fc )| 2 = , §2 _ m 2 )2 T F 2 m 2 4(g 2 - ^5+ - ^L) (43) 

^ ciAy dk dk 

where s,t, u are defined by s = (p e +Puj) 2 = xs, t = (p e — p x +) 2 = —Q 2 , u = rrr x+ +m 2 k — s—t 
and the coupling constants 7, 5, e are 

im dk U 12 



V2M W cos 



im dj U 12 

8l = 1l , 5r = -^=— 1 -. (44) 

V 2M W cos p 



We follow here the notation of f55 |, where one can find the expressions for the matrices 
Uij, Vij which diagonalise the chargino mass matrix. 
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9.3 Virtual Squark t-channel Exchange at LEP 

The differential cross-section can be expressed as 

da{e Tl q " ) -P M i + ^-^ +A ^ < 45 > 

acosu acosu Sins 

and the amplitudes are given by 

A, = X ' 4 t 2 
1 4(t-m 2 ) 2 

Ate Q e Qq „ r 9z a e a q i\ / Aa \ 

M = ~ (T^) +Ue{ s-Ml + iT z M z }) (46) 

Ax, A 2 correspond to the contributions from the t-channel diagram and the SM interference 
respectively. Here fh is the mass of the exchanged squark; Q e , Q q are the electric charge of 
the electron and quark q respectively; gz = sin0 e cosd ; is the angle between the incoming 
electron and the quark q; t = — 1(1 + cos^); and coupling constants defined by 

a e = -^ + sin 2 6 w 
a u = ---sin 6 W 

a d = ^sin 2 ^ (47) 
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